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July - December 1993 Activity under NASA Contract NAG8-974 


FINANCIAL SUMMARY 

Through December 31 , 1993 our spending on this contract totals $14, 379.83 out of 
$100,000 received. We anticipate spending to increase considerably over the remaining 6 
months of the contract. However, we foresee a request for a one-year no-cost extension which 
will allow us to refine our techniques with promising data sources anticipated from GOES-I. 

SCIENTIFIC SUMMARY 

The thrust of the proposed effort under this contract is aimed at improving techniques to 
track water vapor data in sequences of imagery from geostationary satellites. In regards to this 
task, significant testing, evaluation and progress was accomplished during this reporting 
period. Some of the results are compiled in a paper presented at the 2nd International Satellite 
Winds Workshop (Appendix 1). To summarize, sets of winds derived from Meteosat data were 
routinely produced during Atlantic hurricane events in the 1993 season. These wind sets were 
delivered via Internet in real time to the Hurricane Research Division in Miami for their 
evaluation in a track forecast model. For eighteen cases in which 72-h forecasts were 
produced, thirteen resulted in track forecast improvements (some quite significant). In addition, 
quality-controlled Meteosat water vapor winds produced by NESDIS were validated against 
rawinsondes, yielding an 8 m/sec RMS. This figure is comparable to upper-level cloud drift 
wind accuracies. Given the complimentary horizontal coverage in cloud-free areas, we believe 
that water vapor vectors can supplement cloud-drift wind information to provide good full-disk 
coverage of the upper-tropospheric flow. The impact of these winds on numerical analysis and 
forecasts will be tested in the next reporting period (outlined below). 



A new water vapor motion target selection procedure has been developed and tested for 
use in cloud-free areas. The motivation for the development is that water vapor motion tracers 
are not necessarily related to pixel brightness, unlike cloud tracers. The maximum gradient 
(MG) procedure uses a moving window of size X (pixel space) to select a target by finding the 
maximum gradient (bi-dimensional pixel value differencing) within the window. The gradient is 
tested for all combinations of pixels within a pre-defined distance in the target window. Multiple 
target selection is optional. The MG procedure differs from current techniques that select 
targets from clouds by eliminating the search for maximum brightness within the target box. 
This is allowable given the fact that the water vapor channels do not sense land pixels or 
gradients, as is the case in IR window or visible channel tracking. Our evaluation of the new 
method indicates an increase in targets and vectors of around 25%, with a reduced vector 
RMS compared to collocated rawinsondes of .4 m/sec. The MG target selection procedure has 
been incorporated into the CIMSS wind tracking system for 'clear-air' water vapor tracer 
determination. 


The coverage, accuracy and impact of water vapor motion vectors will depend on an 
optimum balance between clearness of image (signal to noise), resolution of sensor, and time 
interval between images. This relationship was investigated using special dwell sound modes 
from the VAS instrument on-board GOES-7. Data (imagery) was obtained from both large 
(16km) and small (8km) detectors. The large detector imagery was derived from a single 
sample (spin), however the small detector imagery was derived from multiple sampling (4- 
spin), thereby reducing considerably the signal-to-noise ratio. Water vapor vectors were 
derived for two spectral bands (6.7 and 7.3 micron), and compared statistically to collocated 
rawinsondes. The greatest impact occurs with the 7.3 micron data, which normally carries a 
high signal-to-noise ratio. The combination of improved horizontal resolution and reduced 
noise results in significant improvement in vector quantity (by 40-50%) and quality (25% 
improvement in vector RMS). Tracking accuracy is less affected by image time separation. It 



was found that 1-h intervals are slightly preferred over 2-h or 30-minute, although 30-min. 
intervals were slightly superior for the 6.7 micron small detector data. Further experiments are 
currently being conducted using GOES-7 water vapor and IR data. 

The final step in the wind tracking algorithm is objective editing. Improvements in the so- 
called CIMSS 'auto-editor' during this reporting period include: 1) tuning of the penalty function 
to minimize RMS errors with rawinsondes, 2) addition of a maximum-allowable height re- 
assignment, 3) incorporation of an empirically-based bias adjustment, and 4) additional 
flexibility added to the code to allow for improved inspection and evaluation. This version of the 
auto-editor is being implemented operationally at NESDIS. We plan on investigating the 
optimization of the objective editor for water vapor winds as part of this contract. 

We are continuing to investigate improved methods for vector height assignments. In 
addition to fine tuning the C02 slicing algorithm (which is currently operational at NESDIS), an 
H20-intercept method is being evaluated in anticipation of GOES-I (which will not contain a 13 
micron C02 channel on the imager). A comparison of the two methods indicates they yield 
similar results (Nieman et al. 1993, JAM, 32, 1559-1568), although the C02 method appears 
slightly superior. Further refinement of both methods is currently being investigated, or is 
planned. 

Collaboration with NASA-MSFC on the issue of optimizing satellite wind tracking algorithms 
was deemed an essential component of our proposed effort. Previous studies have suggested 
that the CIMSS and MSFC algorithms may contain independent strengths and weaknesses, 
and an effort to evaluate these and combine the attributes of each system is desirable. To this 
end, software routines developed at CIMSS were delivered to MSFC for incorporation and 
evaluation. Specifically, routines developed to calculate superior vector height assignments 
using multi-spectral information (C02 data) were passed to MSFC. It is hoped an 



intercomparison of algorithms can be conducted in the future to address the optimum tracking 
algorithm issue. 

CIMSS is also collaborating with NASA-MSFC in the comparison of satellite-derived wind 
vectors with wind observations from jimsphere and radar wind profilers. Preliminary results 
indicate good agreement between the jimsphere and satellite data (4.4 m/sec RMS), and an 
effort is underway to enlarge the sample with the intention of formal publication. 

During the month of March (1994), CIMSS will conduct a data collection exercise in 
collaboration with the European Space Operations Center (ESOC), which will be jointly 
sponsored by the GOES-I project. Water vapor wind sets will be produced twice a day from 
GOES-7 (CIMSS) and Meteosat-5 (ESOC). This exercise will serve several purposes: 1) 
provide a good data base for satellite intercomparison studies, 2) demonstrate the potential for 
the wind sets to benefit projects such as GEWEX and PATHFINDER, 3) provide a benchmark 
for future wind-tracking algorithm testing and evaluation, and 4) serve as a data set for model 
impact studies. NMC and ECMWF have been contacted and have expressed an interest in 
near-real time acquisition and assimilation of the data into their respective global assimilation 
systems. 

PUBLICATIONS 

Velden, C.S., 1993: Investigation of water vapor motion winds. Postprints of 2nd International 
Satellite Winds Workshop, Tokyo, Japan, 13-15 December. (Appendix 1) 

Velden, C.S., and S.J. Nieman, 1994: Investigation of water vapor motion winds from 
geostationary satellites. To be delivered at the 7th Coni. Satellite Meteor., Monterey, CA, 6-10 
June. (Abstract in Appendix 2) 
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ABSTRACT 


Motions deduced in animated water vapor imagery from geostationary satellites can be used to infer wind fields 
in cloudless regimes. For the past several years, CTMSS has been exploring this potentially important source of 
global-scale wind information. Recently, METEOSAT-3 data has become routinely available to both the U.S. 
operational and research community. Compared with the current GOES satellite, the METROS AT has a superior 
resolution (5km vs. 16km) in its water vapor channel. Preliminary work at CIMSS has demonstrated that wind sets 
derived from METEOSAT water vapor imagery can provide important upper- tropospheric wind information in data 
void areas, and can positively impact numerical model guidance in meteorological applications. Specifically, 
hurricane track forecasts can be improved. Currently, we are exploring methods to further improve the derivation 
and quality of the water vapor wind sets. 


1. INTRODUCTION 

Water vapor imagery from GOES satellites has been available for over a decade. These data are used 
extensively, mainly in a qualitative mode, by forecasters in the United States (Weldon and Holmes 1991). Some 
attempts have been made at quantifying the data by tracking features in time sequences of the imagery (Stewart et 
al. 1985; Hayden and Stewart 1987). For a variety of reasons, applications of this approach have produced marginal 
results (Velden 1990). Recently, METEOSAT-3 (M-3) was repositioned at 75W by the European Space Agency, in 
order to provide complete coverage of the Atlantic Ocean. Data from this satellite are being transmitted to the U.S. 
for operational use. Compared with the GOES satellite, the M-3 has a superior resolution and signal- to- noise ratio in 
its water vapor channel, which translates into improved automated tracking capabilities. 

During the Atlantic hurricane seasons of 1992 and 1993, tracking algorithms were applied to the M-3 water 
vapor data in order to evaluate the coverage, accuracy and model impact of the derived vectors. Data sets were 
produced during several tropical cyclone cases, including Hurricane Andrew. In this paper, the M-3 water vapor 
wind sets are assessed, and their impact on a hurricane track forecast model is examined. 

2. STATISTICAL EVALUATION OF THE M-3 WATER VAPOR WINDS 

During the spring of 1992, M-3 water vapor wind sets were routinely produced on a daily basis (around 
1200UTC) at CIMSS. The domain of the wind sets covered the eastern U.S. and the western North Atlantic ocean 
basin (Fig. 1). The purpose of this exercise was to assess the horizontal and vertical coverage of the winds, and 
evaluate their accuracy. M-3 water vapor targets were selectively chosen to be constrained to approximate cloud- 
free areas in the imagery. This was accomplished by activating an empirically-determined brightness temperature 



threshold value which was not to be exceeded in the target selection step of the wind-tracking algorithm. An 
objective quality control procedure is invoked once the wind vector field has been derived (Hayden and Velden 
1991). 

The data set presented in Fig. 1 is typical of the data sets produced during the exercise. From a purely qualitative 
point of view, the horizontal coverage of the water vapor wind vectors shown in Fig. 1 is quite good in comparison 
to conventional and cloud-drift wind observations routinely available over the western North Atlantic basin. It was 
found that the vertical distribution of the assigned vector pressure -heights in cloud-free areas (based on brightness 
temperature to radiosonde temperature comparisons) was typically in the range of 200-500mb, with a peak near 
350mb. It was also clearly demonstrated for future considerations that the water vapor wind sets could be created on 
McIDAS (or VDUC at NMC) in a time scale commensurate with real time operations. 



Fig. 1. Typical M-3 water vapor wind set coverage. 


Comparisons between the M-3 water vapor wind vectors and collocated rawinsondes (within 1.0 degree) over the 
period April - June 1992 were compiled and are presented in Table 1. The winds are also evaluated against the 
collocated first guess forecast (in this case the Aviation model 12h forecast). Both vector speed bias and RMS were 
computed and compared. The data sets are also stratified to reflect comparisons with selected rawinsondes in 
relatively remote areas. 


Table 1. Statistical evaluation of M-3 water vapor motion winds from the period April- June, 1992. 


1) Versus eastern US/westem 
Atlantic rawinsondes (N=981) 


2) Versus Bermuda rawinsondc 
(N=31) 


3) Versus Guadeloupe rawinsonde 
(N=50) 


H20 Winds 

NMC Forecast 

H20 Winds 

NMC Forecast 

H20 Winds 

NMC Forecast 

Speed bias 
(m/s) 

-0.9 

-1.7 

-0.8 

-2.1 

-0.7 

-1.9 

Vector RMS 
(m/s) 

7.8 

6.6 

6.9 

7.1 

6.6 

6.8 


Overall, the M-3 vector speed bias (wind vector minus rawinsonde) is -0.9 m/sec, which is just under 1 m/sec 
better than the first guess forecast. The RMS error, however, is 1.2 m/sec higher than the first guess. The superior 
performance of the first guess RMS error can be explained by the fact that most of the comparisons were over the 
eastern United States, an area where the model has been properly initialized with abundant rawinsonde data. 




Examination of the relatively remote Bermuda and Guadeloupe rawinsonde comparisons, however, show the RMS 
error of the M-3 winds slightly below that of the first guess. The speed bias is also much improved over that of the 
first guess. 

More recently, M-3 wind sets have been routinely produced on a daily basis by NESDIS in Washington DC since 
May 1993. The algorithms and tracking methods were transferred from CIMSS in order to test and evaluate them in 
an operational environment The target selection and thresholds, and quality control procedures are nearly identical 
to those utilized at CIMSS. Validation against collocated rawinsondes (N=1471) indicates a vector RMS of 8.2m/s, 
and a speed bias of -1.1 m/s, agreeing closely with the results in Table 1. 

3. IMPACT ON NUMERICAL HURRICANE TRACK FORECASTS 

Another way of quantitatively evaluating new data types is through model impact studies. In our investigation, a 
hurricane track forecast model (VICBAR, DeMaria et al. 1992) is used to test the sensitivity of the water vapor wind 
data on numerical hurricane track forecasts. VICBAR is a nested, spectral barotropic model that has been run in 
near-real time at the NOAA/Atlantic Oceanographic and Meteorological Laboratory Hurricane Research Division 
(AOML-HRD) for the past few years. The initial condition for the forecast model is a vertically averaged (mass 
weighted) deep- layer mean wind over the 850-200mb depth of the troposphere, with an added tropical cyclone 
bogus. The barotropic forecast model uses the shallow-water equations, with the forecast storm track determined 
from the location of the relative vorticity maximum on the innermost model mesh. This forecasting system has been 
used to evaluate impacts of other data types (Franklin and DeMaria 1992; Velden et al. 1992) 

Wind sets were produced daily at 1200UTC during several Atlantic tropical cyclones in 1992 and 1993. For the 
model impact evaluation, forecasts were included based on the following criteria: 1) the storm was within the 
domain of the generated wind set, and 2) the storm intensity was tropical storm strength or better (>35kts). There 
were 32 forecast cases, from 10 different storms (including Hurricane Andrew in 1992) that met these criteria. 

Table 2 shows the results of forecasts which were initialized with operationally available data plus the water 
vapor wind set produced at CIMSS, compared to control runs which were initialized on operationally available data 
only. Considering all of the cases together, the results show modest improvements to the forecasts with the inclusion 
of the water vapor wind data. The improvement in mean forecast error (MFI in Table 2) ranges from 2.2% at 24h to 
6.8% at 72h. After an adjustment is made for the serial correlation between forecast cases (Franklin and DeMaria 
1992), none of the improvements were found to be significant at the 95% confidence level. Another measure of 
forecast impact is the frequency of improved forecasts (FIF), which simply shows the percentage of forecasts which 
resulted in some improvement when the water vapor data were included. The FIF at 72h indicates that 71% (15 out 
of 21) of the forecasts were improved. 


Table 2. Impact of METEOSAT water vapor winds on VICBAR tropical cyclone track forecasts. The following 
verification statistics are valid for forecasts which included the water vapor winds: number of forecasts 
(N) , effective number of independent forecasts (N*), and the mean forecast error (MFE) relative to best 
track verification. Also given are comparisons with the control forecasts: the mean forecast improvement 
(MFI), expressed in both kilometers and as a percent relative to the control forecast error, the standard 
deviation of the improvements (SOI), the number of improved forecasts (IF), the frequency of improved 
forecasts (FIF) expressed as a percent, and whether the forecast improvements are statistically 
significant at the 95% confidence level (SIG) . 


Forecast 


Interval 

(h) 

N 

N* 

MFE 

(km) 

MFI 

(km) 

MFI 

% 

SDI 

(km) 

IF 

# 

FIF 

% 

SIG 

(y/n) 

24 

32 

26.6 

149.1 

3.3 

2.2 

23.8 

19 

60 

n 

48 

27 

22.4 

338.4 

3.9 

1.1 

48.2 

16 

60 

n 

72 

21 

17.4 

347.0 

23.6 

6.8 

91.6 

15 

71 

n 


Several of the selected tropical cyclone cases were relatively close to the U.S. mainland. It is reasonable to 
assume that in these cases, the VICBAR model was relatively well initialized by the nearby conventional 
(operational) data base, limiting the potential for the satellite data to have a positive impact. On the other hand, 
storms well out to sea should make better candidates for forecast improvement. To test this hypothesis, the sample 
was stratified to include only those cases east of 70W. The results of these 23 cases are shown in Table 3. While 
still not statistically significant at the 95% confidence limit, the MFI percentages and FIF are a notable 
improvement over the statistics presented in Table 2. 


Table 3. Same as Table 2, except only for cases east of 70W. 


Forecast 


Interval 



MFE 

MFI 

MFI 

SDI 

IF 

FIF 

SIG 

(h) 

N 

N* 

(km) 

(km) 

% 

(km) 

# 

% 

(y/n) 

24 

23 

19 

136.4 

9.0 

6.6 

20.2 

16 

70 

n 

48 

20 

17.8 

335.3 

10.2 

4.3 

60.1 

13 

65 

n 

72 

18 

16.2 

362.3 

30.4 

8.4 

98.6 

13 

72 

n 


Using the stratified sample, it is of interest to examine the distribution of forecast differences between the runs 
made with the water vapor winds and the control runs in order to fully appreciate the FIF results. Table 4 shows the 
distribution of the 72h forecast differences. Of the eighteen cases, only two forecasts that included the water vapor 
winds were notably degraded (greater than 100km) relative to the control forecasts. One of these poor forecasts was 
a Hurricane Bonnie case, which seemed to result at least in part because Bonnie’s track did not follow the deep layer 
flow used in the model forecast. Rather, it followed a shallow-layer flow, in a direction quite different from the 
deep-layer flow. The water vapor winds provide measurements in the 200-500mb layer, thus affecting the upper part 
of the deep layer mean wind flow field. In this case, the control analysis approximated the shallow-layer flow more 
closely, and as a consequence, the VICBAR model control forecast without the winds was closer to the observed 
storm track. The reason for the other poor forecast has not yet been identified. 

On a more positive note. Table 4 shows that thirteen of the eighteen 72h forecasts were improvements over the 
control forecasts, with five of them being notable improvements (defined here as greater than 100km). Three of the 
notably improved forecasts occurred during Hurricane Andrew’s interaction with an upper-level cyclonic circulation. 
An example is shown in Fig. 2, from 1200 UT 19 August 1992. At this time the upper-level circulation was situated 
to the north of Andrew and quite evident in the water vapor imagery. The derived water vapor winds captured the 
circulation, as shown in Fig. 1. During the next 72h, Fig. 2 shows that Andrew's track was only slightly affected by 
the upper-level low, shifting it north-northwest for a short time, before Andrew escaped its influence and turned 
more to the west However, the VICBAR control forecast from this time recurved Andrew to the north and 
eventually to the northeast in response to the influence of the upper low on the deep layer mean steering flow. 
Although the turn back towards the west was not predicted, the VICBAR forecast made with water vapor wind data 
responded with much less curvature, and an improved longer-range forecast. 

Table 4. Distribution of 72h forecast errors relative to control forecast errors for cases east of 70W. 

Forecast differences (km) 

* * * * ********** * * * * 

I - I I I -I 

200 100 0 100 200 
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Fig. 2. Example showing impact of water vapor winds on a VICB AR forecast of Hurricane Andrews track. 

4. WORK IN PROGRESS 

To this point, the CIMSS water vapor winds have been produced from automated techniques designed for cloud 
drift wind extraction. This includes image preprocessing, target selection, the tracking algorithm, height assignment, 
and quality control. Clearly, all five of these basic procedures can be refined to better suit the characteristics of the 
water vapor data. Efforts are underway at CIMSS to examine these five areas to better tune them to water vapor 
tracking, and some examples are briefly discussed below. 

a) Image preprocessing 

NASA-MSFC has developed image filtering techniques to help reduce the noise in the water vapor imagery 
before the target selection process is invoked. These techniques will be investigated in collaboration with MSFC. 

b) Target selection 

The current target selection procedure involves identifying the coldest pixel in the target box with an 
accompanying brightness temperature gradient which satisfies certain thresholds. This philosophy was chosen to 
maximize the likelihood of selecting cloud edges as targets, since these features represent good tracers over time in 
conventional IR imagery. Since we are attempting to provide vectors in cloud-free areas when using the water vapor 
data, this philosphy may not be optimal. A method will be investigated that simply searches for the maximum 
brightness temperature gradient over the target selection box. The resultant targets will be evaluated against 
routinely produced targets. 

c) Tracking algorithm 

The current automated tracking algorithm employed at CIMSS (Merrill 1989) uses a correlation method to 
identify potential matches of target radiance features in subsequent imagery. Currently, a sequence of three images, 
each separated by one hour, are used to derive two vectors. The two vectors are then averaged to achieve a final 
vector which is used to represent the tracer. Preliminary investigations indicate the averaging philosphy may not be 
appropriate with water vapor tracking since, compared to cloud tracking, correlations are lower and noise is higher. 
An algorithm will be investigated that selects the vector (first, second or average) that best satisfies accompanying 
quality control checks. In addition, the use of a four image sequence instead of three will be studied. The four image 
sequence will yield three vectors, which would provide a better evaluation of tracer consistency. 

d) Height assignment 

The method currently employed at CIMSS simply matches the target brightness temperature with a first guess 
temperature profile (assumes a one to one relationship) and extracts a pressure height. A statistical evaluation of the 
actual relationship may yield superior vector height assignments (Hayden and Stewart 1987). 



e) Quality control 


A two-step quality control procedure is used in the CIMSS tracking system. The first step involves consistency 
checks between the first and second vector. The second step is an elaborate automated editor which analyzes the 
wind fields in conjunction with the first guess, and outputs confidence weights on each vector. The auto editor also 
attempts to find the optimum pressure height assignment for each vector, and therefore may reassign the original 
height to one which better fits the first guess profile. Efforts are underway to refine this system to take into account 
the water vapor vector characteristics. 

To optimize the extraction of quality motion vectors from METEOSAT water vapor imagery, CIMSS is 
collaborating with NASA-MSFC and EUMETSAT, and will also entertain ideas from the workshop. Model impact 
experiments are planned in 1994 in coordination with NMC, NASA-MSFC, ECMWF and AOML-HRD. 

5. SUMMARY 

The recent availability of METEOSAT-3 data over the western Atlantic Ocean has led to an effort to extract 
quantitative information from the water vapor channel. Data sets containing vectors derived from animated M-3 
water vapor imagery were produced during 1992 and 1993 using the CIMSS automated wind derivation algorithm. 
A statistical evaluation of the vectors reveals that the water vapor winds (relative to collocated rawinsondes) show a 
reduced speed bias compared to the collocated first guess forecast values. The vector RMS errors are larger (by 
about 1 m/sec) than the first guess over the eastern U.S., but become slightly lower than the first guess at remote 
locations (e.g., Bermuda and Guadeloupe). 

The water vapor wind sets were also demonstrated to have a slightly positive impact on barotropic numerical 
hurricane track forecasts (VICBAR model), although the results were not statistically significant at the 95% level. 
Except for a couple of examples, the preliminary results seem to suggest that the positive impact is maximized on 
cases well offshore and away from data-rich regions, as would be expected from intuitive reasoning. Most of the 
cases near the U.S. coast showed negligible or slightly negative impact. The most notable forecast improvements 
occurred during Hurricane Andrew’s interaction with an upper-level circulation, which was well -captured by the M- 
3 water vapor winds. Overall, 72h track forecasts were improved by 7%, while nearly 3/4 of the VICBAR forecasts 
showed some improvement with the inclusion of the water vapor winds into the initial analysis. 

Current and future efforts are being directed towards refining the CIMSS wind extraction and quality control 
procedures to better account for the water vapor data and resultant motion vector characteristics. It is envisioned 
that these data can compliment cloud drift wind data to provide complete global up per- tropospheric wind coverage. 
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Motions deduced in animated water vapor imagery from geostationary 
satellites can be used to infer wind fields in cloudless atmospheric regimes 
For the past several years, CIMSS has been exploring this potentially 
important source of global -scale wind information. Recently, Meteosat data 
has become routinely available to both the U.S. operational and research 
community. Compared to the current GOES satellite, the Meteosat has a 
superior resolution (5km vs 16km) in its water vapor channel. Taking 
advantage of this, preliminary work at CIMSS has demonstrated that wind sets 
derived from Meteosat water vapor imagery can provide important 
upper- tropospheric wind information in data void areas, and can positively 
impact numerical model analyses and guidance. One specific application is in 
hurricane track forecasting. In this paper, the Meteosat water vapor motion 
wind sets are assessed, and their impact on numerical analysis and forecast 
models is examined. Current work at CIMSS continues to explore methods to 
further improve the derivation and quality of the water vapor motion winds, 
and this will be discussed at the conference. 



